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Abstract 



In this thesis we study the prospects for observing the invisibly decaying Higgs boson in the 
associated tiH production at the LHC. The results of the Monte Carlo simulations of signal 
and background processes show that there is a possibility of observing the statistically 
significant number of signal events required for the discovery. Moreover, the analysis can 
be further improved to reduce the number of false reconstructions of the W boson. 

The analysis of the tiH production is independent of the model in which the Higgs boson 
decays into the invisible channel. There are several possibilities for models where H 
invisible can be of interest. For this thesis, we have studied the simplest supersymmetric 
model, called mSUGRA. The results of the scans of the mSUGRA model parameter space 
show that the regions, where the branching ratio of the lightest neutral Higgs boson to 
the lightest neutralino pair is high, are excluded by current experimental constraints. The 
h ^ bb channel dominates, and the possibility for discovery in this channel will not be 
suppressed by the invisible decays. This result does not disqualify invisible channel as 
possible signature in other models. 

Chapter m of this thesis is based on the published paper: B.P. Kersevan, M. Malawski . 
E. Richter-Wqs: Prospects for observing an invisibly decaying Higgs boson in the ttH 
production at the LHC, The European Physical Journal C - Particles and Fields, 2003, vol. 
29, no. 4, pp. 541 - 548. 
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Chapter 1 
Introduction 



The Higgs boson is a very important element of the modern theory of elementary particles 
and their interactions, called the Standard Model. By means of spontaneous symmetry 
breaking the Higgs mechanism is used to introduce mass to other particles. The Higgs 
boson remains the last undiscovered particle predicted by the Standard Model. The most 
important large-scale experiments of these days, at the Large Hadron Collider at CERN 
and at Tevatron at Fermilab aim at the discovery of this particle. 

Worldwide collaborations of scientists are studying the possible production and decay 
channels, where the Higgs boson may be observed. There are theories, in which the Higgs 
boson can decay into the invisible channel. This can occur when the Higgs decays into 
particles which interact very weakly with the matter, so they cannot be observed in the 
detector. There are various models that may lead to such invisible decays. These models 
include the decay into the lightest neutralinos in the supersymmetry models, or into neu- 
trinos in the models of the neutrino mass generation, such as extra dimensions, TeV-scale 
gravity, Majorana models or 4th generation lepton [T2] . 

The associated production of the Higgs boson and a tt-quark pair can be an interesting 
channel leading to the observation of the invisible Higgs boson. The top quark production 
is a process with a very characteristic topology, so it can be used to effectively filter the 
events in which Higgs production takes place. One of the objectives of this thesis was to 
study the possibility of observing the Higgs boson in this channel using simulations the 
Atlas detector at LHC. 

Supersymmetry (SUSY) is considered by many theorists as the most probable exten- 
sion of the Standard Model. SUSY assumes a symmetry between fermions and bosons, 
introducing superpartners for each of the known elementary particles. The Minimal Super- 
symmetric extension of the Standard Model (MSSM) predicts the existence of the Lightest 
Supersymmetric Particle (LSP), which does not decay into normal particles and is a can- 
didate for dark matter. A Higgs boson, decaying into a pair of LSPs, can give an invisible 
signature. 

In order to obtain a more realistic estimation of the tiH process analysis results, it is 
necessary to discern the cross-section of the tiH production and the branching ratio of the 
Higgs to the invisible channel. SUSY models allow for such predictions. In this thesis, we 
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describe the results of the scans of the mSUGRA model parameter space. The objectives 
of the scans were to find such regions in the parameter space, where the branching ratio 
of the Higgs into the invisible channel is large enough to make the results of the analysis 
statistically significant. The production process was studied only for the lightest Higgs 
boson in the so-called decoupling limit, i.e. in the region where its coupling to fermions 
and bosons is almost equal to that of the Standard Model Higgs boson of the same mass. 

The thesis is organized as follows: Chapter [21 provides a theoretical introduction to the 
Standard Model Higgs boson and gives an overview of the current and planned experiments 
for Higgs boson discovery. Chapter [HI provides a brief introduction to supersymmetry and 
the specific mSUGRA model that can lead to an invisible Higgs decay. In Chapter HI we 
describe the analysis of the tiH, H invisible process based on Monte Carlo simulations 
of the ATLAS detector at LHC. ChapterElpresents the results of the scans of the mSUGRA 
model parameter space, where the invisible decay channel is possible. The conclusions are 
summarised in Chapter El 
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Chapter 2 



The Higgs Boson 



2.1 The Standard Model and the Higgs Theory 

The Standard Model (SM for short), is currently regarded as a successful theory describing 
elementary particles and their interactions. SM is based on the quantum field theory. 
It assumes the existence of fermion fields representing matter and gauge fields that carry 
interactions. Matter is composed of three generations of quarks and leptons. For the spinor 
fields ip, we have the Dirac Lagrangian: 

L = ^{i-1^'df,-m)%lj. (2.1) 

>From this Lagrangian we can derive the Dirac equation, which is the equation of motion 
of free relativistic particles and antiparticles (fermions). In order to introduce interactions 
into the model, we postulate the gauge invariance of the fields. The Lagrangian should be 
invariant under the U{1) gauge transformation: 

^ ^ e-^^^'^V- (2.2) 

For the Lagrangian (|2.H) to be invariant, we need to introduce a new vector field and 
derive a new Lagrangian: 

Ld = L- e^YAf^iP, (2.3) 
where the transformation law for the vector field is: 

A^^ A'^ = A^ + ^d^aix). (2.4) 

To the Lagrangian ()2.3|1 we also add a term that yields the equations for gauge fields. For 
the electromagnetic field, the term is: 

La = --F^.F^^\ F,, = d^A, - d,A,. (2.5) 
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Such a Lagrangian can provide the Maxwell equations for electromagnetic fields and the 
Dirac equations for fermions: 

nA'' = 6-07^ = f; l^iid^j, - eA^)'ijj -mijj^O. (2.6) 
In the case of a multiplet of N particles, we have a more general Dirac Lagrangian: 

L = ^{iYdf, - M)^, (2.7) 

where \l/ is an iV— component vector. When we require the gauge invariance of such a 
theory, the transformation has the following form: 

= G*, (2.8) 

where G is the N x N unitary matrix, detG — 1. These matrices form the group SU{N). 
The SU (N) invariance requires the introduction of N'^ — 1 vector fields and modification 
of the Lagrangian: 

Ln = L-g^YT.T''G;^: (2-9) 

n 

where are generators of the SU{N) group. By introducing the operators: 

G'm = EG'm^"' G^.u = d^G,-d,G^ + ig[G^,G,], (2.10) 

n 

we can construct the Lagrangian for the vector fields: 

LG^~Tr{G^,Gn (2.11) 
The full Lagrangian of the SU{N) theory has the following form: 

Lym = ^mid^ - gG^^l - - lTr{G^,Gn (2.12) 

To describe the electroweak interactions for one generation of leptons or quarks, we 
need the SU{2)l x U{1) group. The doublet ^ consists of the (lepton, neutrino) pair or 
(up-, down-) type quark pair respectively. The left-handed components ^z, = L of these 
spinors transform as the fundamental representation of the SU{2) group, while the right- 
handed components = R are singlets with respect to this group. All members of 
the doublet transform under the U{1) group. Consequently, we need four vector fields, 

corresponding to U{1) and Sj, {k = 1, ...,3) for SU{2). These fields give the fermion 
Lagrangian: 

Lf = iLYd^,L-L^^{\g'YLB^ + gB';T^)L 

+ ii^],Yd,il^l-tl^],Y\g'Y^RB,,lj], (2.13) 
+ t^lrd.i^l-i^lrlg'Y^RB.ijl. 

In this Lagrangian, the weak hypercharge is denoted by Y. The Lagrangian for vector 
fields has the following form: 
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Chapter 2. The Higgs Boson 



(2.14) 



The fermion Lagrangian (I2.13|] cannot contain the mass term which is not invariant 
with respect to the SU{2) group, since the left and right components transform differently, 
and: 



= ^^^^ + 



(2.15) 



Moreover, fields are also massless, since we do not obtain the linear mass terms in 
motion equations from the Lagrangian ()2.14|) .. 

The model with all particles being massless is unsatisfactory, since we observe massive 
particles in the experiments. The solution to this problem is the Higgs mechanism. We 
add to our model one more SU{2) doublet, consisting of two complex scalar fields: 



H 



hi 
h2 



(2.16) 



called the Higgs field. The Lagrangian for these fields has to be gauge-invariant under U{1) 
and SU{2). We can choose such a gauge in which the Higgs field has the form: 



H 




h 



(2.17) 



where h is the real function. This means that there is only one physical degree of freedom, 
so the Higgs field can give one neutral scalar particle. Very important in the construction 
of the Higgs Lagrangian is the substitution of the mass term by the Higgs potential. We 
postulate: 



(2.18) 



When fi^ < the Higgs potential has a non-trivial minimum as shown in Fig. 12.11 The 
minimum is for: 



H^H -- 

We can choose the phase of the solution as: 



H = Ho 





v/V2 



2X' 



2 



(2.19) 



(2.20) 



Such a solution breaks the gauge symmetry. This is called the spontaneous symmetry 
breaking, since the potential itself is gauge-invariant, whereas the vacuum state breaks the 
symmetry. The Lagrangian containing interactions of the Higgs and the vector fields has 
the following form: 
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Figure 2.1: Mexican-hat Higgs potential with a non-zero minimum 



^ k 



(2.21) 



This Lagrangian contains the terms bilinear in fields. We need to diagonalize this expres- 
sion and derive mass coefficients. For components Bj^ and B"^ the Lagrangian is already 
diagonal, so we obtain the mass terms: 



= = — 
2 

We can form combinations for particles with a definite charge: 

= — 

These states form the charged W bosons that carry the electroweak interactions. We also 
have combinations that yield the neutral bosons: 



f2.22) 



(2.23) 



-g'B, + gBl 



(2.24) 



The mass of these bosons is: 



The remaining combination diagonalizing ()2.2H1 is 



A.. 



g'B.+gBl 



(2.25) 



(2.26) 



This component is not present in the Lagrangian, so the mass coefficient must vanish. The 
massless field represents the photon. The ratio between the masses of W and Z bosons 
defines the Weinberg angle: 
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w 



cos 9w- 



(2.27) 



Mz Vg''' + g' 

To construct the Lagrangian for the Higgs field, we can rewrite it by distinguishing the 
vacuum part from the dynamical part: 



H{x) ^ Ho + h{x), h{x) 





h{x) 



(2.28) 



The Lagrangian for this field can be written in the form: 



Lh = d^,hd^h + 2ji^h'^ -2\/2Xvh^ - Xh^ + —^[vh + h\Z^Z^, + 2cosewW^W-]. (2.29) 

COSnw ^ ^ 

This Lagrangian gives the coupling of the Higgs field to the vector bosons and the Higgs 
mass: 



m,H — —2fj,^ — 2Xv. 



(2.30) 



This means that the value of the v parameter is not enough to determine the Higgs mass, 
but it also depends on A, which is a free parameter of the Standard Model. 

The Higgs mechanism also gives mass to fermions. This is done by introducing the 
Yukawa couplings to the Lagrangian: 



where 



^ = [/lo + h{x)] 



1 




For h{x) = we obtain 



(2.31) 



(2.32) 



(2.33) 



giving us the fermion masses: 



mi = Giho] m,2 = (^2^0- (2.34) 
Similarly, we obtain the couplings of the Higgs field to fermions in the Lagrangian: 



-Gih{x)ipi'il)i - G2h{x)'ijj2ilJ2- 



(2.35) 



>From these equations we can see that the coupling of the Higgs boson to fermions is 
proportional to the fermion mass. As a consequence, we can predict that the Higgs boson 
will mostly decay into the heaviest lepton or quark pairs. 
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2.2 Current Experimental Results 

The Standard Model predictions have passed all experimental tests conducted so far. The 
electroweak gauge bosons W and Z were detected at the LEP experiment at CERN and 
the t quark was observed at the Tevatron. Only the Higgs boson remains as the last 
undiscovered elementary particle of the Standard Model. Discovery of the Higgs boson 
will enable verification of theoretical foundations of the Standard Model. 

The possibilities of observing the Higgs particle both in e+e~ and in hadron colliders 
have been extensively analyzed. Since the coupling of the Higgs to other particles is 
proportional to their masses, the Higgs production and decay often involve the heavy 
quarks t or b. Sample processes of Higgs production, such as Higgs-strahlung and weak 
boson fusion are shown in Fig l2.2[ 




Figure 2.2: Processes of Higgs production 



The branching ratios of the Higgs boson to the others particles depend on the Higgs 
mass. Fig. 12.31 shows the branching ratios for various possible decay channels. 
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Figure 2.3: The branching ratios of the Standard Model Higgs boson 

The results gathered at LEP from direct search for the Higgs boson and from the 
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precision electroweak measurements provide some approximate limits to the Higgs boson 
mass. Direct search at LEP has excluded the existence of the Higgs boson with a mass 
rriH < 114.4 GeV, at a 95% confidence level. The Fig l2.4l shows the expected and observed 
behavior of the statistic defined in [l]. 



^ 50 
^40 
^30 
20 
10 

-10 



-20 



106 108 110 112 114 116 118 120 

mH(GeV/c^) 

Figure 2.4: The results of direct search for the SM Higgs boson at LEP [T]. 

The precise measurements of the electroweak boson masses and couplings can also give 
an estimation of the Higgs boson mass. The combined results from the LEP and Tevatron 
experiments are shown in Fig. 12.51 They suggest the preferred Higgs boson mass to be: 
rriH = 113 (+62 — 42) GeV at a 68% confidence level. They also set the upper limit to be 
about 237 GeV at a 95% confidence level. 

2.3 LHC Discovery Potential 

The LHC experiment will search for the SM Higgs boson in collisions of proton-proton 
pairs. The ATLAS and CMS experiments are preparing for possible discovery channels. 
Fig. 12.61 shows the ATLAS discovery potential for the SM Higgs boson for the integrated 
luminosity of 30 fb~l |2]. It can be seen that the total estimated significance is higher than 
the 5(7 discovery threshold for the Higgs mass range between 100 and 200 GeV. 




Observed 

Expected for background 
Expected for signal plus background 

J \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 
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Figure 2.5: The fit to electroweak measurements at LEP, setting bounds on the Higgs 
boson mass [T6] 
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Figure 2.6: ATLAS discovery potential for the SM Higgs boson |21 
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Supersymmetry 



3.1 Introduction to Supersymmetric Theories 

In spite of their success and compatibility with experimental results, both the Standard 
Model and the Higgs mechanism outlined in Chapter [21 remain unsatisfactory as theo- 
ries [lOj. One of the problems that the Standard Model does not solve is the hierarchy 
problem. This is a problem of the sensitivity of the Higgs mass to quantum corrections 
from the particles coupling to the Higgs boson. It so happens that the corrections from 
bosons and from fermions have opposite signs. Thus, the corrections would cancel each 
other out if each of the known fermions had a bosonic partner of the same mass and each 
boson had a fermionic partner. Such a symmetry is called supersymmetry |19| . 

Supersymmetry transforms boson states into fermions, and vice versa. This can be 
expressed as follows: 

Q\Boson) = \Fermion); Q\Fermion) = \Boson) (3-1) 

The supersymmetry generator Q is an anticommuting spinor that satisfies the following 
algebra: 



{Q,Qt} = PM 
{Q,Q} = {Q\Q^} 



[P", Q] 









(3.2) 
(3.3) 
(3.4) 



Supersymmetric particles can form multiplets, consisting of bosonic and fermionic states, 
called supermultiplets. All particles and their superpartners from the same supermultiplets 
have the same masses, since the operator commutes with the Q, . The supersymmetry 
generators commute also with the generators of gauge transformations. Hence, particles in 
the same supermultiplet have the same electric charge, weak isospin and color degrees of 
freedom. In each supermultiplet, the number of fermionic degrees of freedom np is equal 
to the number of bosonic degrees of freedom ub- 



11 



Section 3.2. Minimal Supersymmetric Standard Model 



Up = ub (3.5) 

The simplest possibility for a supermultiplet includes a single Weyl fermion with two 
helicity states and a complex scalar field with two degrees of freedom. Such a combi- 
nation is called the chiral or matter supermultiplet. The next simplest possibility is the 
supermultiplet with a massless spin-1 vector boson, which has two degrees of freedom. Its 
superpartner is a massless spin-1/2 Weyl fermion. Such a combination is called the vector 
supermultiplet. 

3.2 Minimal Supersymmetric Standard Model 

The Standard Model particles and their superpartners are presented in Table 13.11 The 
table shows the states before the electroweak symmetry is broken, i.e. the W and B bosons 
correspond to and B^ fields from Lagrangian ()2.14jl and their partners are called winos 
and binos. After breaking this symmetry, new states appear: zino Z and photino 7. All 
these particles and their interactions form the Minimal Supersymmetric extension to the 
Standard Model (MSSM). 

In Table 13.11 we can see two Higgs doublets: H^, H^. In the Standard Model, one 
complex scalar Higgs doublet was enough (see eq. (|2.16|) ). Now, however, two complex 
doublets are required to give mass to the particles. 

Table 3.1: The Standard Model particles and their superpartners before electroweak sym- 
metry breaking 



Particle 




Spin 


Spartner 




spin 


quarks: 


u, d 


1/2 


squarks: 


u, d 





lepton: 


I 


1/2 


sleptons: 


I 





Higgs: 


Hu, Hd 





higgsino: 


Hu, Hd 


1/2 


gluon: 


9 





gluino: 


9 


1/2 


W bosons: 


W 





winos: 


W 


1/2 


B bosons: 


B 





binos: 


B 


1/2 



The free Lagrangian of the supersymmetric model can be written in the following form: 

%ree = -9''(t>*'d^<P, - t^^'a^dfi^i + F*'F,, (3.6) 

where are the scalar fields, ifj are the Weyl spinors from the supermultiplet and F is the 
auxiliary scalar field, required to fill the remaining scalar degrees of freedom. 
The interaction can be introduced into the theory by adding the Lagrangian 
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knt = -l^^i^j + W'F, + c.c, (3.7) 

where W''^ and W' are the functions of the bosonic fields (pi. These terms can be expressed 
using the superpotential W: 

W = \m''Mj + \y''^M,(i>k, (3.8) 

and relations: 

ocpiocpj d(pi 

where M*-^ is the mass matrix for the fermion fields, and y''^^ is a Yukawa coupling between 
the two fermions ipii^j and the scalar 0^. 

In the case of the MSSM, the superpotential has the form: 

^MSSM = ^VuQHu - dyaQHd - eyeLHd + fxHuHd. (3.10) 

In this equation we use the superfield notation, i.e. Hu, Hd, u, etc. denote the chiral 
supermultiplets, as in Tab. 13. II The term corresponds to the Higgs mass in the Standard 
Model. 



3.3 Models of Supersymmetry Breaking 

The experimental results show that supersymmetry has to be broken, since no sparticles 
with masses equal to their Standard Model partners have been discovered. Despite the 
symmetry breaking, we would like to preserve the supersymmetric relations between su- 
perpartners in the high energy sector of the theory, in order to keep the solution to the 
hierarchy problem valid. Such an approach is called the soft supersymmetry breaking. 
This can be expressed by an effective Lagrangian of the MSSM: 

L^55M _ -l{M3gg + M2WW + MiBB) + c.c. 

- {uauQHu - dadQHd - eaeLHd) + c.c. 1 1 \ 

~ ~ ~ ~ ~+ - -t ~ ~+ W.iij 

- Q'^m'^Q — L^m\L — um\u — dmjd — em\e 

- ml^HlHu-m\H*dHd-hHuHd + c.c.). 

We can see that in eq. (j3.11ll we have the Mi, M2 and M3 gaugino mass terms, a^, a^, ae are 
the matrices of Yukawa couplings, and the mf terms are mass matrices of fermion families. 
Together with the Higgs sector parameters, it is necessary to have 105 free parameters in 
the MSSM. The experimental results that restrict the fiavor-changing neutral currents and 
CP violation suggest that the squark and slepton mass matrices should be diagonal. This 
organizing principle helps reduce the number of parameters. 
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Such universality can be achieved by assuming that, at high energy scale, there exists 
some simple Lagrangian in which there are only a few common mass and coupling param- 
eters. After applying the renormalization group equations to the masses and couplings 
to compute their values in the electroweak scale, we can obtain the full spectrum of the 
MSSM parameters. Such unification of couplings is in agreement with the GUT theories, 
which predict unification at the mass scale of lO^^GeV. Since the high-energy scale theory 
cannot be observed directly, it is called the hidden sector, while the MSSM is called the 
visible sector. 

There is more than one model that can lead to such supersymmetry breaking. All these 
models can be divided with respect to the type of interaction between the hidden sector and 
the visible one. In gravity-mediated sypersymmetry breaking one assumes the supergravity 
Lagrangian in the hidden sector and the gravitational interactions lead to supersymmetry 
breaking. On the other hand, in the gauge-mediated supersymmetry breaking models, one 
assumes that there exist some additional chiral supermultiplets, called messengers, which 
couple to the fields of the hidden sector and also to the MSSM fermions through gauge 
interactions. 

An especially interesting model, called mSUGRA, is the minimal version of gravity- 
mediated supersymmetry breaking models. It assumes unifications of all fermion masses 
into mi/2, all boson masses into mo, and one common coupling constant Aq. To complete 
the model, it is necessary to add the tan/? and sign /i parameters describing the Higgs 
sector. Thus, there are only 5 parameters: 



that can give the full MSSM spectrum after applying renormalization group equations to 
the electroweak scale. 

The mass spectrum of the MSSM is subject to the electroweak symmetry breaking, 
similar to that of the Standard Model, which leads to mixing of the states. In the Higgs 
sector, we have two complex scalar doublets = {H:^,H^) and Hd = {H^,H^). The 
Higgs potential has a minimum for non-zero = and = v^, and the ratio between 
them is called tan /3: 



When the electroweak symmetry is broken, the Higgs doublets mix with the gauge bosons, 
to give masses to the and Z° bosons. Thus, from the eight original degrees of freedom 
only five remain. They form the five scalar boson states: CP-odd neutral scalar A*^, two 
charged scalars H'^,H~, and two CP-even scalars /i° and It can be shown that h*^ is 
the lightest among these Higgs bosons of the MSSM. 

Electroweak symmetry breaking causes the mixing among the sparticles and produces 
new mass eigenstates. Neutral higgsinos H^, and neutral gauginos B, produce four 
neutral states called neutralinos denoted x° In a similar way, the charged higgsinos and 
winos form four charged states, called charginos - two with positive and two with negative 
electric charges. 



mo, mi/2, Ao, tan /3, sign^, 



(3.12) 



tan/3 = Vu/vd 



(3.13) 
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In the MSSM we can define a new multiplicative quantum number, called i?— parity. It 
is defined as: 

= (_l)3{B-L)+2s^ ^3_;L4) 

where B is the barion number, L is the lepton number, and s is the spin of the particle. The 
Standard Model particles and the Higgs bosons have even i?— parity (Pr = 1), while all 
sparticles have odd i?— parity (Pr = —1). The conservation of the i?— parity has important 
phenomenological consequences. It implies that the lightest supersymmetric particle (LSP) 
with Pr = — 1, must be stable, and all sparticles must eventually decay into a state which 
contains an odd number of LSPs. 

One of the best candidates for the LSP is the lightest neutralino Xi- Since it is electri- 
cally neutral and interacts only weakly with ordinary matter, it becomes a good candidate 
for dark matter required by cosmology. Also, in the collider experiments, particles decaying 
into neutral LSPs can give an invisible signature (missing energy) in the detector. 

3.4 Experimental Signals for Supersymmetry 

Although supersymmetry helps solve many theoretical and phenomenological problems 
beyond the Standard Model, there is still no experimental evidence that supersymmetry 
really exists. No sparticles have been discovered so far. The LEP and Tevatron experiments 
have only set limits to the masses of supersymmetric particles [TT]. E.g. the lightest Higgs 
boson has been excluded in the mass range below 91.0 GeV. 

There are plans to search for supersymmetry at the LHC. Most of the sparticles with 
masses below ITeV may be within reach of this collider. Fig. 13.11 shows the regions in 
the (m^,tan/3) MSSM parameter space, where the various Higgs bosons can be observed. 
Most probable is the detection of the neutral h boson, which is the lightest of the Higgs 
particles. Fig. 13.21 shows the regions, where h can be observed in the tih, h ^ bb channel. 
The possibility to observe the lightest Higgs boson in the invisible decay is studied in the 
following chapters of this thesis. 
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Figure 3.1: LHC discovery potential for the MSSM Higgs bosons [TT 
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Figure 3.2: The Sa-discovery contour curves for the tt and Wh with h ^ bb channel in the 
{rriAytanP) plane for the ATLAS experiment at luminosities of 30 fb~^ 0. 
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Chapter 4 

Analysis of UH Production 



4.1 Introduction 

In this chapter we describe the analysis of the process that may allow us to observe the 
Higgs boson in the invisible channel. We describe the specific process and study possible 
backgrounds. In subsequent sections, we will provide details on the selection criteria that 
were applied to extract the interesting signal events from the background. The results and 
possible improvements are presented at the end of this chapter. 

The possibilities to observe the Higgs boson in the invisible channel have been dis- 
cussed throughout the past years. The prospects of detecting such a Higgs particle via its 
associated production with a gluon, Z ov W bosons are described in jn|. The exploitation 
of associated ti production was suggested in [12]. The weak boson fusion process leading 
to the production of invisible Higgs is discussed in IH]. A more recent evaluation of the 
production with gauge bosons can be found in [T2] . 

In this thesis we study the associated tiH production as suggested in [13]. The proposed 
approach was to require the leptonic decay of one of the W bosons and the hadronic decay 
of the second one, as shown in Fig. 14. 11 The final products are a single lepton, 2 6-jets and 
2 jets from the hadronic decay of the W boson. These jets and leptons can be identified 
by the detector so they can be used to trigger the filtering of interesting events. The 
evidence of the invisibly-decaying Higgs boson will be an excess of the selected events over 
the predicted background. 

4.2 Signal and Background Processes 

In order to be able to observe such an excess of signal events, it is necessary to simulate and 
analyze the signal and possible backgrounds. The signal was simulated using the PYTHIA 
[2T] Monte Carlo event generator. The generator simulated the proton-proton collisions 
at the LHC, at the 14 TeV center-of-mass energy, 
gg, qq ttH : 

This is the signal process generated using the PYTHIA event generator. The simulation 
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Figure 4.1: Process for observing the invisible Higgs. 



is configured in such a way that 100% of the Higgs particles are decaying into the invisible 
channel, meaning they are not observed by the detector. The simulation assumes the 
Standard Model coupling of the Higgs to the t quarks. 

The background may come from processes which give a signature that is similar to the 
one generated by the signal process, i.e. an isolated lepton, 2 identified 6-jets, 2 additional 
jets and missing transverse energy. The processes that may lead to such signatures, are 
the t-quark or 6-quark production associated with the W- or Z-boson. 
gg, qq tt: 

This is the main background process. It can yield the same signature as the signal, 
but the main difference is in the missing energy. In the case of the signal, it contains the 
energy of the invisible Higgs, while in the background it mostly carries the energy of the 
neutrino from the leptonic decay of the W boson. Consequently, the missing energy will 
be the used to reduce this background. 

gg, qq ttz, z^uv : 

This background can yield the same signature as the signal. The neutrinos coming from 
the Z decay can carry significant missing energy, comparable with the signal, 
qq ^ ttW, W ^ ^1/ : 

In this process, there can be 3 bosons: one from the associated production of the 
tt pair and 2 from the top-quarks decays. As a consequence, the results of the simulation 
where only the 1^-boson from the associated production is forced to decay into a lepton 
and the neutrino, should be finally multiplied by a combinatorial factor of 3. It assumes 
that the acceptance of the analysis is independent of the source of the lepton. 
gg, qq ^ bbZ/7*, Z/7* ^ © jets : 

This process can yield the same signature as the signal and its cross-section is several 
orders of magnitude larger. However, it is possible to reduce this background considerably 
by requiring the reconstruction of the top-quark from the hadronic jets. Also, requiring 
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significant missing energy and vetoing additional leptons will suppress this background, 
qq bbW, W — > © jets : 

This background can be suppressed by requiring the reconstruction of the top quark in 
the hadronic mode and significant transverse energy. 

4.3 Event Generation and Detector Simulation 

For the tt process event generation, the PYTHIA [21] and HERWIG HH] Monte Carlo 
generators were used. Additionally, the AcerMC [IH] event generator was used to simulate 
other backgrounds. For these generators the initial parameters were set to default values. 

For detector simulation, the ATLFAST [213] program was used. It implements a sim- 
plified simulation of the ATLAS detector at the LHC. The generated events serve as input 
to this program and it produces the detector response as output. It provides such observ- 
ables as the 4-momenta of isolated leptons, jets, tagged 6-jets and transverse energy. It is 
also possible to evaluate the missing transverse energy by integrating the total energy of 
particles detected and subtracting from the initial energy of colliding protons. 

The simulation assumes 90% efficiency in lepton identification, 80% for jet reconstruc- 
tion and 60% for 6-jet identification (tagging). These efficiencies were taken into account 
while calculating the normalized number of events. 

4.4 Analysis 

The objective of the analysis is to apply such filters to the simulated processes that will 
suppress, as much as possible, the background while preserving a high number of signal 
events. The following filters are applied to the events: 

1. One isolated lepton is required. The lepton comes from semi-leptonic decay of one 
top quark. Events with more leptons are rejected in order to eliminate events coming 
from the bbZ/j* background. 

2. Two identified (tagged) 6-jets are required. 

3. One top quark is required to be reconstructed from the hadronic decay mode t — > jjb. 
First, only such jj pairs are selected from all possible combinations of jets, where 
rrijj = rriwi^ 15 GeV. Additionally, only central jets with pseudo-rapidity < 2.0 
are allowed to reduce the number of events with jets coming from initial- or final-state 
QCD radiation instead of the VT-boson decay. Next, the masses of all combinations 
of the selected jj pairs and 6-jets are calculated. The four-momenta of the jets 
are recalibrated so that the rrijj gives rriw exactly. The best combination of the jjb 
system is chosen and the top quark is considered reconstructed iirrijjh = mt±25 GeV 
(see Fig. 112)). 
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Figure 4.2: Reconstructed mass of the top quark for the ttH events (number of events not 
normalized). 



4. It is not possible to reconstruct the second top-quark which decays in the leptonic 
channel. This would require precise information on the missing energy, which in the 
case of the tiH events comes both from the leptonic decay of the top-quark and from 
the invisible Higgs. Instead, it is possible to use the transverse mass of the lepton 
end system to distinguish between the signal and the ti background, where 

mrii, ^t) = \j2E^TpT{l-cos(t)) (4.1) 

The Fig. 14. 31 shows the distribution of mT{i, ^t) for signal and background processes. 
It can be seen that for the ti process, where the missing energy comes from the 
leptonic decay of the VT-boson, there is a sharp end of the distribution around the W 
mass. To cut the big fraction of the ti events we require that mT{i, ^t) > 120 GeV. 

5. In addition to the mT(£, ^t) cut, we also require a large missing transverse energy 
of the event: fir > 150 GeV. 

6. The signal-to-background ratio is enhanced by the additional requirement of large 
transverse momenta in the reconstructed system, X)Pt"^ > 250 GeV. The YjVt^ = 
Y.Pt^Vt where the sum runs over the transverse momenta of reconstructed objects 
from top-quark decays: two 6-jets, two light jets used for the reconstruction of the 
W ^ qq decay and an isolated lepton. This further suppresses the background where 
true top quarks are not present, like hbZ and hbW . 

7. Finally, further enhancement of the signal-to-background ratio can be achieved by 
the additional requirement regarding cone separation, Rjj, between jets used for the 
W jj reconstruction, the Rjj < 2.2. 
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i?,, = ^{Ar^y + (A0)2 (4.2) 
is a distance in the {rj, cj)) plane between the jets. 




Figure 4.3: Reconstructed transverse mass of the lepton and system in the ttH events (left 
plot) and in the events (right plot). The dashed line denotes the distributions calculated from 
the true invisible energy of the primary products of W boson decays in these events, obtained by 
using the generator level information. The distributions are normalized to the number of events 
expected for an integrated luminosity of 30/6~^. 



Table 14.11 shows the cumulative acceptance of the specified filters applied to selected 
processes. 

We can see that the acceptance of signal events is about 0.3%, while the background is 
reduced by a factor of 10^5 — 10~6. 

After performing such selections, one can see that about 70% of the ti events come 
from the lepton-tau decay and 20% from the lepton-lepton decay of the top-quark pair in 
the PYTHIA sample with compatible fractions also found in HERWIG events. The lepton-tau 
label denotes one top quark decaying t — > Wb iub and another t — > Wb rub, where 
i stands for electron or muon. The lepton-lepton decay labels events with both top quarks 
decaying t — > Wb — > iub. Finally, the lepton-hadron decay labels events with one top quark 
decaying t — > Wb — > iub and another t — > Wb qqb. 

In the cases of lepton-lepton and lepton-tau decays, the jjb combination is thus made 
not from the true W ^ qq decays, but from other jets coming from the initial- or final-state 
radiation (ISR/FSR). These events could hopefully be suppressed further by implementing 
a tau-jet veto and with more stringent requirements in the t jjb reconstruction. The 
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Table 4.1: The cumulative acceptances for the specified selection criteria. Efficiencies for fa- 
tagging and lepton identification are included. A Higgs boson mass of 120 GeV is assumed for 
signal events. Only the dominant background sources are listed. 



Process 


ttri 


4-+ V 

ttZj 


it 


tt 






PYTHIA 


AcerMC 


PYTHIA 


HERWIG 




Trigger lepton 


22% 


22% 


22% 


22% 




2 b-jets + 2 jets 


5.0% 


4.8% 


4.9% 


5.2% 




rec. t-quark (jjb) 


2.6% 


2.4% 


2.4% 


2.6% 




m^^^^ > 120 GeV 


0.87% 


0.93% 


4.1 • 10"^ 


5.2 • 10" 


-4 


> 150 GeV 


0.41% 


0.53% 


2.3- 10-5 


3.7-10- 


-5 


J2pT > 250 GeV 


0.40% 


0.51% 


2.0- 10-5 


3.2 • 10" 


-5 


< 2.2 


0.28% 


0.35% 


7.5 • 10-6 


1.2 • 10" 


-5 



signal events, in contrast, contain only a ~ 10% fraction of lepton-tau and lepton-lepton 
decays. The relative fractions of signal and ti background events are shown in Figure 

Considering the relative fractions of the tau-lepton events in the signal and background, 
one can assume that the inter-jet cone separation, might provide some additional 
separation power; the for signal and tt background events are given in Figure 14.41 
Subsequently, a loose cut of < 2.2 was applied; the final efficiencies are listed in 
Table lO 

The expected numbers of events for an integrated luminosity of 30/6-^ are given in 
Table 14.21 Several values of the Higgs boson masses are studied, while assuming the 
Standard Model production cross-sections. The branching ratio of the Higgs boson to the 
invisible channel is assumed to be 100% in this analysis. 

It can be seen that for the process tt generated with PYTHIA, the signal-to-background 
ratio is about 39% for a Higgs boson mass of 100 GeV and about 9% for a Higgs mass of 
200 GeV. Of course, if only true reconstructed top-quark decays were allowed, these ratios 
would increase considerably. 

It would also be possible to increase the signal-to-background ratio by setting more 
tight cuts on such observables as -Rjj, J2Pt^- However, such optimizations are not 
suitable at this moment, since they could lead to over-fitting to the results obtained by 
the particular implementations of the processes in event generators. It is clear that further 
optimization of the analysis should lead to an increase in the efficiency of rejecting "fake" 

~^ jj reconstructions and eliminating the contribution of events with lepton-tau and 
lepton-lepton decays. 
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Figure 4.4: The relative fractions of the tt decay modes are hsted for signal and tt background 
simulated with PYTHIA (left plot). The cone separation between jets used in the W jj 
reconstruction; the distributions are normalized to the number of events expected for an integrated 
luminosity of 30 fb~^ (right plot). 



The basic measure of the quality of the analysis is the significance S = s/Vb^ where s 
and b are the numbers of signal and background events respectively. 5 = 5 is assumed to 
be the discovery level, while = 3 is called the evidence level. S = 1.96 corresponds to 
the 95% confidence level. We can define the parameter as: 

e^ = -^xBr(H^inv). 

It can be used to estimate the minimal branching ratio required for the significance to 
reach the desired thresholds. 

Taking the values from Table IT!2l we get s = 60 for the ttH signal with tjih = lOOGel^ 
and 6 = 65 for the total background that includes the (lep-had) ti events generated with 
PYTHIA. This yields S = 7.44 and requires to be 0.67 and 0.4 for the discovery and 
evidence levels respectively. The values for the ti background and also for (lep-had) events 
are shown in Table 
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Table 4.2: Expected numbers of events for an integrated luminosity of 30fb~^ and selection as 
specified in Table ITTl Efficiencies for b-tagging and lepton identification are included. The (PY) 
and (HW) denote the results for the ti events generated with PYTHIA and HERWIG respectively. Also 
shown is the separate contribution to the tt background from the lepton-hadron events. 



Process 


No. of events 


tiH, 




ruH = 100 GeV 


60 


mjy = 120 GeV 


45 


mjy = 140 GeV 


30 


mj? = 160 GeV 


25 


mH = 2m GeV 


15 


tiZ 


20 


ttw 


20 


tt (all) 


115 (PY) , 190 (HW) 


(only lepton-hadron) 


15 (PY) , 30 (HW) 


bbW 


5 


bbZ/j* 


5 



Table 4.3: Estimated significance 5* = s/Vb for tt background and estimated branching 
ratios required for significance to reach the values of 5, 3 and 1.96 





all tt 








lep-had tt 








ruH 


s/Vb, 






^2(1.96) 


s/Vb 








100 GeV 


4.67 


1.07 


0.64 


0.42 


7.44 


0.67 


0.4 


0.26 


120 GeV 


3.5 


1.43 


0.86 


0.56 


5.58 


0.9 


0.54 


0.35 


140 GeV 


2.34 


2.14 


1.28 


0.84 


3.72 


1.34 


0.81 


0.53 


160 GeV 


1.95 


2.57 


1.54 


1.01 


3.1 


1.61 


0.97 


0.63 


200 GeV 


1.17 


4.28 


2.57 


1.68 


1.86 


2.69 


1.61 


1.05 
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5.1 Objective of the scans 

The supersymmetric models as described in ChapterElallow a decay of the Higgs boson into 
a pair of the lightest neutralinos XiXv This is an invisible channel when the neutralino is 
the LSP, which means it does not decay, but leaves the detector, leaving its signature only 
as missing energy. It is interesting to study for which parameters of the supersymmetric 
models, like mSUGRA, such decay is allowed and can lead to possible Higgs discovery. 
Results of such scans can be combined with results from similar scans for Higgs decaying 
into visible Standard Model particles. 

In this thesis, the Suspect and HDECAY |Hl[Z| programs were used to perform the scan of 
the mSUGRA model parameter space. The goal of the scan was to identify the regions with 
high Higgs branching ratio to XiXi pair. The analysis of the associated tiH production 
presented in |0] and discussed in Chapter |3| was performed with the assumption that the 
BR{H — i> XiXi) was 100%. Therefore, it is necessary to multiply the results of the analysis 
by the branching ratio obtained from the scans. The estimates for minimal branching ratios 
required for the discovery or evidence signatures respectively are included in Table In 
the scans, we will search for such areas in the parameter space, where the BR{H XiXi) 
can reach these values. 

It is also interesting to check the branching ratio of Higgs to the bb pair. This decay is 
considered as the main channel where the Higgs can be observed in the same production 
mode. If the decay of if — > XiXi is open, the H ^ bb will become suppressed. We can 
compare the branching ratio to the bb and the invisible channels and examine if the invisible 
decay should be considered as complementary for the observation. 

5.2 Implementation of the Scans 

For the scans, we have combined two programs. Suspect [8j takes as input the 5 parameters 
of the mSUGRA model (see eq. 1)3.12^ land returns the spectrum of masses of MSSM 
particles. These parameters serve as input to the HDECAY [Zj program, which computes the 
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branching ratios of the Higgs to both Standard Model and MSSM particles. 

The input parameters were selected to cover the most interesting regions of the mSUGRA 
space. The masses rriQ and mi/2 were covered in the range of to 500 GeV, with steps 
of 10 GeV. The tan (3 had values of 50, 25 and 12.5, . . . Three values of the Aq parameter 
were chosen: -200, 0, 200 and both values of sign /i were used. 

As the input for the Suspect and HDECAY programs, default parameters were used. For 
the drawing of plots, the ROOT framework was used |3j. 

5.3 Constraints on the Parameter Space 

It is important to note that not all regions of the mSUGRA parameter space are allowed. 
First, it is possible that the mass of the stau is lower than that of the lightest neutralino. 
This would imply that the neutralino is not the LSP, contradicting the assumption that 
SUSY particles decay into the invisible channel. Such regions should be theoretically 
excluded. 

Second, current experimental results from the LEP and Tevatron experiments, as well 
as from cosmological observations have excluded the existence of SUSY particles below 
certain mass limits. Such limits for the SUSY Higgs boson obtained by four combined 
LEP experiments is now set at 91.0 GeV [l^. More limits on the masses of other sparticles 
can be found in [TT]. The constraints applied to the scans are shown in Table l5Tl 

Additionally, the programs that were used for the scans have limitations related to 
numerical stability. Although several tests are performed internally by these programs, 
there are also points of the parameter space, where the calculations do not converge. Such 
points are also excluded from the scans. 

Table 5.1: Experimental constraints on masses of sparticles used in the scan 



Particle 


Mass constraint 


h 


> 91.0 GeV 


sleptons 


> 80.0 GeV 


chargino 


> 103.6 GeV 


gluino 


> 195.0 GeV 


neutralino 


> 45.6 GeV 



5.4 Results of the Scans 

The resulting plots showing the contours of branching ratios are presented in Figs. 15. 1115. 21 
We selected the following values of tan/3: 25, 12.5, 6.25, 3.12; both signs of and produced 
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plots in the (mo, mi/2) plane. We show only values of Aq = 0, since other tested values 
(-200, 200) provide similar results. For tan/3 > 25 and < 3 we observed more instability 
in the combined programs. 

For each value of tan /3 and sign fi we show three contour plots in the {mo, 1711/2) plane: 
branching ratio of lightest Higgs boson to neutralinos, to the bb pair and the Higgs mass 
rrih. The theoretical and experimental constraints are indicated by the solid line: regions 
below and to the left of this line are excluded. 

We can see that there are small regions where the BR{h XiXi) is dominant. For 
the /i > and for small tan/3, the value of the BR can reach 0.9 for small regions of the 
parameter space. This would mean that assuming the possibility of analysis as described 
in Chapter 131a sufficient excess of signal events would be observed to enable observation of 
the Higgs in the invisible channel. However, taking into account the experimental limits, 
it is very likely that the regions where the BR{h XiX\) is large, are excluded by current 
observations. On the other hand, the plots show that the branching ratio to 6-quark is 
dominant in most of the regions that are not excluded. That confirms the need for the 
search for Higgs in this decay channel. 

Figs. 15.31 - 15.51 show the projection of the mSUGRA model onto the (m^, tan/3) plane. 
In the mSUGRA model the niA is not given explicitly as an independent variable, but is a 
function of the initial parameters. It means that one point from the {niA, tanjd) plane may 
correspond to more than one (^0,^1/2) pair of initial parameters. In such a case, when 
plotting the values of branching ratios, we chose the maximum from these values. 

Fig. 15.31 shows the branching ratio BR{h XiXi) iii ^ mSUGRA model without any 
theoretical and experimental constraints. We can see that there are regions where the 
XiXi channel can be dominating. However, when we apply the constraints according to 
Table Em we can see in Fig. 15.41 that all points with a high branching ratio are excluded, 
which is also seen in the (mo, mi/2) plane. On the other hand, the BR{h — > bb) can be 
very high even with constraints applied (Fig. I5.5|l . 
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Figure 5.1: Scans of mSUGRA parameter space for yU < 0. 
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Section 5.4. Results of the Scans 



BR(h^XX) H>Ol I BR(h^xx) >i<0l 







Figure 5.5: The maximum values of BR{h bb) in mSUGRA models projected on the 
{mA,tanl3) plane, Aq = 0, after applying theoretical and experimental constraints. 
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Conclusions 



In this thesis we studied the prospects for observing the invisibly decaying Higgs boson 
in the associated tiH production at the LHC. The results of Monte Carlo simulations of 
the signal and background processes have shown that there is a possibility of observing a 
statistically-significant number of signal events required for the discovery. Moreover, the 
analysis can still be improved to reduce the number of false reconstructions of the W boson. 

It is important to note that the analysis of the tiH production described in Chapter HI is 
independent of the model in which the Higgs boson decays into the invisible channel. There 
are several possibilities for models where H — ^ invisible can be of interest. These models 
include the decay into lightest neutralinos in the supersymmetry models, or decays into 
neutrinos in the various models of the neutrino mass generation, such as extra dimensions, 
TeV scale gravity, Majorana models or 4th generation lepton. 

For this thesis we have studied only the simplest supersymmetric model, called mSUGRA. 
The results of the scans of the mSUGRA model parameter space show that the regions 
where the branching ratio of the Higgs particle to the lightest neutralino pair is high, are 
excluded by current experimental constraints. The h —>■ bb channel dominates and the 
prospects for discovery in this channel will not be suppressed by opening of invisible decay. 
This is very encouraging conclusion: we will much prefer to have Higgs boson discovered 
in the visible channel. 

In less constrained SUSY models this however might not be the case. It would be 
also interesting to investigate other models that can lead to such a signature, not just 
the supersymmetric ones. This is out of the scope of this thesis, but might provide an 
interesting research subject. 
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